##### Exerccios 4 - Anlises Exploratrias ###############################################################3333


##### 4.2 Cervejas ##############################################

##Uma amostra de 30 estudantes foi indagada sobre seu tipo de cerveja preferida,com o seguinte resultado 

cervejas <-c("chope","lata","garrafa","chope","garrafa", "garrafa","lata","lata","nenhuma","lata","garrafa","garrafa", "garrafa","lata","lata","lata","garrafa","lata","chope","nenhuma", "garrafa","garrafa","garrafa","chope","garrafa","garrafa","chope","garrafa","lata","lata")
class(cervejas)

cervejas<-factor(cervejas)
cervejas


##Represente este resultado como um grfico de barras e um dotplot (funo dotchart).

par(mfrow=c(2,2))
barplot(table(cervejas))
dotchart(table(cervejas))
par(mfrow=c(1,1))

##Qual tem maior razo dado/tinta?
#Resposta# O grfico dotplot tem a maior razo, pois gasta menos tinta.



##### 4.3 Caixetais ########################################

#Neste exerccio, use o objeto caixeta, criado no tutorial Explorao de uma Varivel Categrica.

caixeta <- read.csv("caixeta.csv", as.is=T)
head(caixeta)
names(caixeta)


#Construa um histograma do DAP dos fustes dos caixetais.

caixeta$dap=caixeta$cap/3.14
caixeta$dap
hist(caixeta$dap, main = "Caixetais DAP", xlab="DAP", ylab = "Frequencia")

#Construa histogramas da altura das rvores para os diferentes caixetais ('local').

library(lattice)
histogram (fuste ~ h | local, data=caixeta, main="Altura Caixetais", xlab="altura", ylab="Frequencia")


#H diferenas entre as estruturas (distribuio de tamanhos) dos caixetais?
#Resposta#  Sim, na Jureia h maior quantidade de indivduos altos e menor de indivduos baixos em relao a Retiro e Chauas.




##### 4.4 Eucaliptos ################################

##Neste exerccio, use o conjunto de dados Inventrio em Florestas Plantadas de Eucalyptus grandis.

egrandis <- read.csv("egrandis.csv", sep = ";", as.is=T)
egrandis
head(egrandis)

##Utilize o grfico boxplot para analisar o DAP de rvores de E. grandis em funo das variveis regio (regiao) e rotao (rotacao).

library(lattice)

class(egrandis$rotacao)
class(egrandis$regiao)

egrandis$rotacao=factor(egrandis$rotacao)

table(egrandis$rotacao)
table(egrandis$regiao)

bwplot(dap~regiao|rotacao, data=egrandis)
bwplot(dap~regiao|rotacao, data=egrandis, main="DAP Eucalyptus grandis por Rotao (1 e 2) por Regio")



##Avalie a normalidade da altura do conjunto total de rvores com um grfico quantil-quantil contra a distribuio normal.

qqnorm(egrandis$ht); qqline(egrandis$ht)
class(egrandis$ht)





##### 4.5 Mais Caixetais#################################

caixeta <- read.csv("caixeta.csv", as.is=T)
head(caixeta)
names(caixeta)


#Aqui usaremos novamente o objeto caixeta, criado no tutorial Explorao de uma Varivel Categrica.

## 1# Analise a relao dap-altura ('dap' e 'h') em funo do caixetal (local) 

#Primeiro fazemos um data frame para "Tabebuia cassinoides"
caixeta.Tab<- subset(caixeta, especie == "Tabebuia cassinoides") 
head(caixeta.Tab) 

#Segundo, calculamos a rea basal
caixeta.Tab$AB=pi*(((caixeta.Tab$cap/3.14)/2)^2)

#Terceiro, somamos os fustes e criamos um dataframe por arvore e por local e por altura
caixeta.arv <- aggregate(caixeta.Tab$AB, by=list(altura=caixeta.Tab$h, arvore=caixeta.Tab$arvore, local=caixeta.Tab$local), FUN=sum)
head(caixeta.arv)

#Quarto, encontramos o DAP atravs da AB
caixeta.arv$DAP=2*sqrt(caixeta.arv$x/pi)
head(caixeta.arv)

#Quinto, relao DAP x altura
coplot(DAP~altura|local, data=caixeta.arv, xlab="Altura", ylab="DAP") 


## 2# Para a mesma relao do item anterior, verifique linearidade com a funo scatter.smooth


caixeta.arv.chauas=caixeta.arv[caixeta.arv$local=="chauas", ]
caixeta.arv.retiro=caixeta.arv[caixeta.arv$local=="retiro", ]
caixeta.arv.jureia=caixeta.arv[caixeta.arv$local=="jureia", ]

# Grficos

par(mfrow=c(1,3))
scatter.smooth(caixeta.arv.chauas$DAP, caixeta.arv.chauas$altura, main = "Linearidade Chaua", xlab = "DAP", ylab="Altura")
scatter.smooth(caixeta.arv.retiro$DAP, caixeta.arv.retiro$altura, main = "Linearidade Retiro", xlab = "DAP", ylab="Altura")
scatter.smooth(caixeta.arv.jureia$DAP, caixeta.arv.jureia$altura, main = "Linearidade Juria", xlab = "DAP", ylab="Altura")
par(mfrow=c(1,1))


## 3# Utilizando o pacote lattice, analise a relao dap-altura ('dap' e 'h') em funo do caixetal (local), mas somente para as rvores de caixeta (Tabebuia cassinoides).

library(lattice)

xyplot (DAP ~ altura | local, data=caixeta.arv, main="Relao Altura x DAP", xlab="altura", ylab="DAP")

